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Energy release due to antineutrino untrapping and diquark 
condensation in hot quark star evolution 
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Abstract. We study the consequences of antineutrino trapping in hot quark matter for quark star configurations 
with possible diquark condensation. Due to the conditions of charge neutrality and /?- equilibrium the flavor 
asymmetry increases with the number density of trapped antineutrinos. Above a critical value of the antineutrino 
chemical potential of 30 MeV diquark condensation is inhibited at low densities and a two-phase structure emerges: 
a superconducting quark matter core surrounded by a shell of normal quark matter. When the quark star cools 
down below a temperature T ~ 1 MeV, the mean free path of antineutrinos becomes larger than the thickness 
of the normal quark matter shell so that they get untrapped within a sudden process. By comparing the masses 
of configurations with the same baryon number we estimate that the release of energy due to the antineutrino 
untrapping transition can be in the range of lO^'^ -r- 10^^ erg. 
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1. Introduction 

The engine which drives supernova explosions and gamma 
ray bursts being among the most energetic pheno mena in 
the universe remains still puzzling [Piran fc Nakail lE)021]. 
The phase transition to a quark matter phase may be 
a me chanism that coul d relea s e such an amoun t of en - 
ergy |Drago fc Tambinil l)l999j) . iBerezhiani et all l)2002|) ]. 
It has been proposed that due to the Cooper instability in 
dense Fermi gases cold dense quark matter shall be in the 
color super conducting st a te with a nonvanishing diquar k 
condensate |Alford et all l|2001^ ■ iBlaschke et all l)200 1a>)]. 
The consequences of diquark condensation for the cooling 
of compact stars due to changes in the transport proper- 
ties and neutri no emissivities have be e n investigate d much 
in detail, see Blaschke et all tOOd^^. IPaee et alJ l|2boQl . 
iBlaschke et al.l l|2001bl) . lBlaschke et all l|2003a|) ]~and may 
even contribute to the explanation of the relative low tem- 
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f)erature of the pulsar in the supernova remmant 3C58 
Grigoria n et all |i002)] . 

Unlike the case of normal (electronic) superconductors, 
the pairing energy gap in quark matter is of the order 
of the Fermi energy so that diquark condensation gives 
considerable contributions to the equation of state (EoS) 
of the order of (A//i)^. Therefore, it has been suggested 
that there might be scenarios which identify the unknown 
source of the energy of 10^^ erg with a release of bind- 
ing energy due to Cooper pairing of quarks in the core of 
a cooling protoneutron star |Hong et al.l l|200ll[) ]. In that 
work the total diquark condensation energy released in 
a bounce of the core is estimated as (A//i)^Mcore corre- 
sponding to a few percent of a solar mass, that is 10^^ 
erg. In this estimate, general rela tivistic effects have bee n 
disregarded. It has been shown in [Blaschke et al.l l|2003a|) ] 
by solving the selfconsistent problem of the star configu- 
rations, however, that these effects due to the stiffening of 
the EoS in the diquark condensation transition lead to an 
increase in the gravitational mass of the star contrary to 
the naive estimates. 
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It has also been demonstrated Blaschke et al.l (|20033)] 
that the energy release due to cooling of a quark core in 
a protoneutron star does not occur whithin an explosive 
process, since the diquark condensation is a second order 
phase transition. 

In the present work, we propose a new mechanism of 
energy release which involves a first order phase transi- 
tion induced by antineutrino untrapping. (Anti-)neutrino 
trapping occurs in hot compact star configurations at 
temperatures T > 1 MeV where the mean free path of 
(anti-)neutrino s is sm aller than the typical size of a star 
" Prakash et al.l l|200ll) and references therein] . 



During the collapse in the hot era of protoneutron 
star evolution, antineutrinos are produced due to the /3- 
processes. Since they have a small mean free path, they 
cannot escape and the asymmetry in the system is in- 
creased. This entails that the formation of the diquark 
condensate is shifted to higher densities or even inhibited 
depending on the fraction of trapped antineutrinos. 

As the quark star cools, a two-phase structure will oc- 
cur. Despite the asymmetry, the interior of the quark star 
(because of its large density) could consist of color super- 
conducting quark matter, whereas in the more dilute outer 
shell, diquark condensation cannot occur and quark mat- 
ter remains in the normal state, opaque to antineutrinos 
for r > 1 MeV. When in the continued cooling process 
the antineutrino mean free path increases above the size 
of this normal matter shell, an outburst of neutrinos oc- 
curs and gives rise to an energy release of the order of 
10^^ — 10^^ erg. This untrapping transition is of first order 
and could lead to an explosive phenomenon. 

The scenario to be detailed in the present paper sug- 
gests that the first pulse of neutrinos emitted in the delep- 
tonization stage of the core collapse, after a cooling time 
scale, shall be followed by a second pulse of antineutrinos 
as an observable characteristics of the present scenario. 



2. Hot 2SC quark matter 

The investigation of the phase structure of electrically and 
color neutral quark matter in /3-equilibrium within quark 
models which a ssume approximate 5't/(3) f lavor symme- 
try ha s revealed Alford fc R.aia gopal 1* 20021) . ISteiner et alJ 
1)2002!) ] that the Color Flavor Locking (CFL) phase is en- 
ergetically favored over the 2-Flavor Superconductivity 
(2SC) one. On th e other hand, i t has been shown 
Gocke et all l|200j) . iNeumann et all l)2003(l ] that at low 
temperatures a sequential deconfinement of light and 



strange quark flavors occurs in dynamical quark models 
which solve the three-flavor gap equations. Thus, strange 
quark matter phases like the CFL one do appear only at 
high densities in the very inner core of compact stars and 
do not occupy a large enough volume in order to cause 
observable effects. Therefore, we consider in the present 
work two flavor quark matter in the 2SC phase only. 



2.1. Thermodynamic potential for asymmetric quark 
matter 

We consider the grand canonical thermodynamic poten- 
tial for 2SC quark matter with in a nonlocal chiral quark 
model [Blaschke et"all ()2003al) ] where in the mean field 
approximation the mass gap (f) / and the diquark gap A ap- 
pear as order parameters and a decomposition into color 
(c G {r,b,g}) and flavor (/ e {u,d}) degrees of freedom 
can be made. 



(1) 



where T is the temperature and /i fc the chemical potential 
for the quark with flavor / and color c. 

The contribution of quarks with given color c and fla- 
vor / to the thermodynamic potential is 



l)^(0/,A;M/c,T)-Fr!S, 



A2 



24 Gi 24 G2 

J / dqq^{oj[ec{Ef{q)+p^fc),T] + 

Jo 

uj[ec{Ef{q) - Hfc),T]} , 



(2) 



where Gi and G2 are coupling constants in the scalar me- 
son and diquark chanels, respectively. The dispersion re- 
lation for unpaired quarks with dynamical mass function 
mf{q) = nif + g{q)(j)f is given by 



Ef{q) = -^q^ + mjiq) . 
In Eq. ||2J) we have introduced the notation 

uj [ec, r] = Tin 



cxp(-- 



where the first argument is given by 



(3) 



(4) 



(5) 



When we choose the green and blue colors to be paired 
and the red ones to remain unpaired, we have 



Ac ^ giq)A{ScM + ScJ 



(6) 



For a homogeneuos system, the thermodynamic poten- 
tial Qq corresponds to the pressure; therefore the constant 
^vac = ^c^vac '^^ choscn such that the pressure of the 
physical vacuum vanishes. 

The nonlocality of the interaction between the quarks 
is implemented via formfactor functions g{q) in the mo- 
mentum space. We use the Gaussian formfactor defined 
as 



g{q) = exp(-gVA2) , 



(7) 



since we coul d show that for t h is choic e stable hybrid stars 
are possible Blaschke et al.l ()2003bf) ]. The parameters: 



A = 1.025 GeV, Gi = 3.761 A^ and m„ = ma 



2.41 



MeV are fixed by the pion mass = 140 MeV, pion de- 
cay constant /^r — 93 MeV and the constituent quark mass 
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mo = 330 MeV at T = = |Schmidt et all (|l994)]. The 

constant G2 is a free parameter of the approach which we 
fix as G2 = 0.86 d. 

We introduce new variables: the quark chemical poten- 
tial for the color c, /igc — (l^uc + fJ'dc)/'^, and the chemical 
potential of the isospin asymmetry, fii = (/i„c ~ fJ-dc)/'2, 
which is color independent. We consider now symmetric 
ifJ-uc = fJ-dc) or nearly symmetric {fi^c — fJ-dc) quark mat- 
ter as the preferable situation to form a pair of fermions 
in the color space, i.e. \^i\ < fiqc- 

The color asymmetry induces a splitting of the quark 
chemical potentials relative to the mean value for two fla- 
vors fiq which is proportional to the new chemical poten- 
tial Us- Therefore we can write. 



5c,g - 2(5c 







(8) 



where /ig and are conjugate to the quark number den- 
sity and the color charge density, respectively. 

We perfom the approximation qi,, = 4)^ = (t> and ex- 



press the thermodynamic potential as in [Kirivama et al.l 

Hooll 



1 f°° 

2 1 dqq^{uj [fLr{-^J'qr " M/): 7"] + 

LU [er{lJ-qr - ^^I),T] + LO [er{-Hqr + ^J.I), T] + 
LO [er{^lqr + ^.l),T]} 

2 

2 I dqq {iv [eb{E{q) - ^,f,) - /i/, T] + 

^ Jo 

uj [ebiE{q) + iiqb) - M/, T] + 
LU [eb{E{q) - iiqb) + T] + 
uj[eb{E{q) + ^.qb) + ^iI,T]} , 



(9) 



where the factor 2 in the last integral comes from the 
degeneracy of the blue and green colors (cb = fg). 

The conditions for the local extrema of fig, correspond 
to coupled gap equations for the two order parameters 4> 
and A 



i=0o,A=Ao 



dA 



= . (10) 

/>=0o,A=Ao 

The global minumum of Clq represents the state of thermo- 
dynamical equilibrium from which all equations of state 
can be obtained by derivation. 



charge neutrality: |7i„ 
also be written as 



= 0, which could 



TLB + nj ~ 2ne — 



color neutrality: rtg = f ^^qr 



^Uqb = 



Here the number densities nj are defined in relation to the 
corresponding chemical potentials fj,j as 

1 

0o,Ao;T 



^^J.j 



(11) 



where the index j denotes the particle species. The baryon 
chemical potential defined as — Sfiq — fij. The solution 
of the color neutrality condition shows that fig is about 5- 
7 MeV in the region of relevant densities (/iq ~ 300 — 500 
MeV). Since the isospin asymmetry is independent of /ig 
we consider fiqc ~ /ig in our following calculations. 

The total thermodynamic potential Q contains besides 
the quark contribution ilq also that of the leptons fi*'^ 



A;//g,^/,/ie,MPe;T') 



qy^, A; fj,q, fll,T) 

J2 f^''(w,T), (12) 



which is assumed to be a massless, ideal Fermi gas 



f]'^(^,T) 



1 



1 
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(13) 
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2.2. Beta equilibrium with electrons and trapped 
(anti)neutrinos 

The stellar matter in the quark core of compact stars con- 
sists of u and d quarks, electrons e and antineutrinos 9^ 
under the conditions of 



/^-equilibrium: d < > u + e 

chemical potentials reads 

Me + Mp. = 



which in terms of 



-2/i/ 



Fig. 1. Mass gap (j), diquark gap A and isospin chemical 
potential /i/ as a function of the baryon chemical potential 
fi B for different values of the antineutrino chemical poten- 
tial fijj^ . Solutions obey /3-equilibrium and charge neutral- 
ity conditions. 

the order parameters </) and A as a function of the baryon 
chemical potential hb at T = 0. Also the solution for the 
isospin chemical potential ni is plotted. We can see that 
with increasing antineutrino chemical potential the 
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absolute value of the chemical potential for isospin asym- 
metry increases (see the inset of Fig.^. We will refer 
to this value as the asymmetry in the system, and when 
it increases the onset of the superconducting phase transi- 
tion is shifted to higher densities. For the region of chem- 
ical potentials fiq where the diquark condensate is nonva- 
nishing its value coincides with that for the case /i/ = 0. 
This region, however, shrinks for /z™*" < . The min- 
imum value of antineutrino chemical potential for which 
the onset of diquark condensation does not coincide with 
the chiral transition is /i™"' = 30 MeV for a density of 
Ais = 1 GeV. 



2.3. Equations of state 

The equations of state are obtained from the thermody- 
namic potential 0. For homogeneous systems, the pres- 
sure is P = — fi, the entropy density is s = —{dQ./dT)^^ 
and the energy density is given by the Gibbs fundamental 
relation 



e = sT - P + HqUq + fljUj. 



(14) 



In Fig. 121 we show the pressure as a function of the baryon 
chemical potential and of the energy density for different 
values of the antineutrino chemical potential. The effect 
of the presence of trapped antineutrinos is twofold: (i) the 
onset of superconductivity in quark matter is shifted to 
higher energy densities with increasing fip^ and (ii) the 
equation of state becomes harder. 
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Fig. 2. Pressure vs. baryon chemical potential (left panel) 
and energy density (right panel) for different values of the 
antineutrino chemical potential fip^ . Due to antineutrino 
trapping the onset of superconductivity in quark matter 
is shifted to higher energy densities and equation of state 
becomes harder. 



3. Hot quark stars and neutrino untrapping 

H ot quark star c onfigurations have been considered first 



m 
in 



Kettner et aI]JT99a] for isothermal configurations and 



Blaschke et al 



for adiabatic ones. When deter- 
mining the mass defect of a compact star configuration 
due to diquark condensation i n a coo ling process, it has 
been shown in 



Blaschke et alJ 1 2003a|) ] that general rela- 



tivistic effects become important. Effects which have been 
estim ated by extr a polati ng the binding energy per Cooper 
pair (Hong et aL 1 ll20nih ] may even be conterbalanced 
when the mass distribution is determined selfconsistenly. 
This observation underlines the necessity to estimate the 
release of binding energy in evolutionary processes of com- 
pact stars which change the equation of state selfconsis- 
tently with corresponding changes in the mass distribu- 
tion an d the gravitational field. In Rcf. [Blaschke et alJ 
however, no trapped (anti)neutrinos have been 
considered and it has been demonstrated that in this case 
the diquark condensation occurs within a second order 
phase transition, so that no sudden release of energy had 
to be expected. The situation may change in the case when 
the transition is inhibited by the presence of trapped an- 
tineutrinos which increase the asymmetry above a critical 
value and therefore prevent diquark condensation. 

3.1. Structure and stability of compact stars including 
trapped (anti)neutrinos 

We consider star configurations which are defined as solu- 
tions of the Tolman-Oppenheimer-Volkoff equations 



dP{r) 
dr 

m(r) 



[e(r) + P{r)][m{r) + inr^Pir)] 
r[r — 2m(r)] 

47r [ r'^e{r')dr' , 
Jo 



(15) 
(16) 



where r is the distance from the center and m{r) the mass 
enclosed in a sphere of radius r. 

The equations are solved for the set of central quark 
number densities rig for which the stars are stable. The 
total mass M = m{R) of the star is defined by the radius 
R being the distance to the star surface which fulfills the 
condition of vanishing pressure P{R) = 0. 

The total mass as a function of radius and central 
quark number density of the configurations is plotted in 
Fig. 13 for the case without neutrinos /ip^ — (sohd hue) 
and for two cases with finite antineutrino chemical poten- 
tial: fj.p^ = 100 MeV (dashed line) and fip^ = 150 MeV 
(dash-dotted line). The latter value relates to our esti- 
mate for the maximum antineutrino number density which 
should occur just at the deconfinement phase transit ion, 
but st ill larger values could also be possible Stcincr e t alJ 
tOO^ ]. From Fig. El we can see that the EoS without an- 
tineutrinos is softer than that with antineutrinos (it has 
a lower pressure at a given energy density) and in Fig|31 
we show that this allows more compact configurations (left 
panel) . Consequentely, the presence of antineutrinos tends 
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Fig. 3. Quark star configurations for different antineu- 
trino cliemical potentials /ip^ — 0, 100, 150 MeV. The 
total mass M in solar masses (Afgun = Af© in the text) 
is shown as a function of the radius R (left panel) and 
of the central number density Uq in units of the nuclear 
saturation density jiq (right panel). Asterisks denote two 
different sets of configurations (A,B,f) and (A',B',f') with 
a fixed total baryon number of the set. 



to increase the mass of the star for a given central density, 
see the right panel of FigO 

In order to estimate the effect of antineutrinos on the 
mass of star configurations we choose a reference config- 
uration without antineutrinos with the mass of a typi- 
cal neutron star Mf — 1.4 Mq, see Fig. O The corre- 
sponding radius is Rj = 10.47 km and the central den- 
sity Uq = 9.69 riQ, where jiq — 0.16 fm^'^ is the satura- 
tion density of nuclear matter. The configurations with 
trapped antineutrinos and nonvanishing /.tp^ to compare 
with are chosen to have the same total baryon number 
as the reference star: Nb = 1-51 Nq^ where Nq is the 
total baryon number of the sun. For = 100 MeV we 
obtain Ma = 1.4025 Mq whereas for 150 MeV 

we have Mb = 1-41 Mq. The differences in the radii are 
Ra- Rf = 0.01 km and Rb - Rf = 0.06 km and the 
changes in the central densities — — —0.016 no 
and Tiq — nf = —0.064 uq, respectively. This is a conse- 
quence of the hardening of the EoS due to the presence of 
antineutrinos. A second set of configurations {A',B',f') 
with a fixed baryon number and Mf = 1.6 Mq is shown 
in the Fig|3 

The mass defect AM^ / ~ Mi — Mf can be interpreted 
as an energy release if there is a process which relates 
the configurations with Mi and Mf being the initial and 
final states, respectively. In the following Subsection we 
discuss a possible evolution scenario and we calculate the 
corresponding energy release. 



3.2. Estimate of energy release in antineutrino 
untrapping 

The antineutrino chemical potential /Ltp^ increases the 
asymmetry of the system and can prevent the diquark 
condensation provided the critical value /Xp^ is exceeded, 
see FigH The inhibition of diquark condensation could 
provide conditions for the explosive release of an amount 
of energy coming from the mass defect introduced in the 
previous subsection. 

During the collapse of a protoneutron star, neutron- 
ization takes place via the inverse /3-process releasing neu- 
trinos which can escape through the yet dilute and cold 
outer shell of the protoneutron star. In the continuation of 
the collapse the temperature and the density increase in 
the protoneutron star core where the proton fraction rises 
again. Antineutrinos are created by the direct /3-process in 
hot and dense matter and cannot escape since their mean 
free path is much smaller than the core radius. When un- 
der these conditions the deconfinement phase transition 
to quark matter occurs, this has two important conse- 
quences: (i) due to high temperatures of the order of 40 
MeV the critical density of the phase transition is much 
lower than at zero temperature and thus a larger frac- 
tion of the star will be in the new phase and (ii) due to 
antineutrino trapping the cooling is delayed. 

The star cools down by surface emission of photons 
and antineutrinos. The region of the star where the tem- 
perature falls below the density dependent critical value 
for diquark condensation, will transform to the color 
superconducting state which is almost transparent to 
(anti)ncutrinos. But nevertheless due to the trapped an- 
tineutrinos there is a dilute normal quark matter shell 
which prevents neutrino escape from the superconducting 
bulk of the star. The criterion for the neutrino untrap- 
ping transition is to cool the star below a temperature of 
about 1 MeV when the mean free path of neutrinos be- 
comes larger than the shell radius Prakash ct al. (200^]. 
If at this temperature the antineutrino chemical potential 
is still large then the neutrinos can escape in a sudden 
outburst. If it is small then there will be only a grad- 
ual increase in the luminosity. An estimate for the possi- 
ble release of energy within the outburst scenario can be 
given via the mass defect defined in the previous subsec- 
tion between an initial configuration with trapped neu- 
trinos (state A or B) and a final configuration without 
neutrinos (state /). In Fig. ^ the integral characteristics 
of the star configuration: the total mass M and the total 
baryon number Nb are plotted as a function of the central 
number density riq. We take two possible initial states (A, 
B with Ma, Mb respectively) with different antineutrino 
chemical potentials and we make the evolution to a final 
state (without antineutrinos, with Mf) at constant total 
baryon number Nb as we show in Fig. |31 The configura- 
tions with trapped antineutrinos (/ip^ = 100, 150 MeV) 
correspond to possible initial states before the outburst 
and the configuration without antineutrinos (/ip^ — 0) 
is the final state after it. The mass defect ranges from 
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Fig. 4. Quark star configurations with diquark condensa- 
tion as a function of the central number density Uq in units 
of the nuclear number density uq . The mass defect for the 
transition from initial configurations with fip^ — 100, 150 
MeV to a final configuration with /ip^ = at constant 
total baryon number is shown. 



AMAf = 0.25 % Mq to AMsf = 1 % M© for the initial 
states A and B, respectively. We see that the release of 
energy can be as big as about 1.5 % of the initial mass. 
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Fig. 5. Mass defect AM and corresponding energy release 
AE due to antineutrino untrapping as a function of the 
mass of the final state Mf. The shaded region is defined 
by the estimates for the upper and lower limits of the an- 
tineutrino chemical potential in the initial state /Zp^ = 150 
MeV (dashed-dotted fine) and /ip^ = 100 MeV (dashed 
line), respectively. 

We can make the same construction for different final 
states and in Fig.|Slwe show the mass defect and the cor- 
responding energy release as a function of the mass of the 



final configuration. Two lines for different initial states 
are plotted: the dashed one for A configurations (with 
/ip^ = 100 MeV) and the dash-dotted one for B config- 
urations (with /ip^ = 150 MeV). The final state / with 
vanishing /ip^ has the mass Mj which defines the con- 
served baryon number involved in the untrapping tran- 
sition. The dots represent the end points of the stable 
configurations. The shaded region in between both curves 
shows the possible mass defect (or corresponding energy 
release) according to our estimate of the range of the an- 
tineutrino chemical potentials considered. 

Since the estimated energy release is as large as 10^^ — 
10^^ erg the scenario suggested in the present paper could 
be discussed as a possible engine driving supernova explo- 
sions and gamma ray bursts. 

A second pulse of antineutrinos is expected to be a 
signal of the untrapping transition. 

4. Conclusions 

We have investigated the effects of trapped antineutrinos 
on the asymmetry and diquark condensates in a quark star 
configurations. By comparing configurations with fixed 
baryon number the release of energy in an antineutrino 
untrapping transition is estimated to be of the order of 
2952 gj-g_ Such a transition is of first order so that an- 
tineutrinos can be released in a sudden process (burst). 
This scenario could play an important role to solve the 
problem of the engine of supernova explosions and gamma 
ray bursts. A second antineutrino pulse is suggested as an 
observable characteristics of the present scenario. 
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